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We study electron correlation effects on quantum criticalities of Lifshitz transitions at zero 
temperature, using the mean-field theory based on a preexisting symmetry-broken order, in 
two-dimensional systems. In the presence of interactions, Lifshitz transitions may become dis- 
continuous in contrast to the continuous transition in the original proposal by Lifshitz for 
noninteracting systems. We focus on the quantum criticality at the endpoint of discontinuous 
Lifshitz transitions, which we call the marginal guantum critical point. It shows remarkable 
criticalities arising from its nature as a topological transition. At the point, for the canonical 
ensemble, the susceptibility of the order parameter x i s found to diverge as lnl/|<5A| when 
the "neck" of the Fermi surface collapses at the van Hove singularity. More remarkably, it 
diverges as |<5A| -1 when the electron/hole pocket of the Fermi surface vanishes. Here 5 A is 
the amplitude of the mean field measured from the Lifshitz critical point. On the other hand, 
for the grand canonical ensemble, the discontinuous transitions appear as the electronic phase 
separation and the endpoint of the phase separation is the marginal quantum critical point. 
Especially, when a pocket of the Fermi surface vanishes, the uniform charge compressibility 
k diverges as \5n\~ , instead of \i where Sn is the electron density measured from the criti- 
cal point. Accordingly, Lifshitz transition induces large fluctuations represented by diverging 
X and/or «. Such fluctuations must be involved in physics of competing orders and influence 
diversity of strong correlation effects. 

KEYWORDS: Lifshitz transition, Fermi surface topology, metamagnetism, charge compressibility, nonan- 
alytic expansion 



1. Introduction 

Lifshitz transition is one of the typical continuous 
quantum phase transitions 1 and has been extensively 
studied for several decades. 2-7 However, despite the in- 
creasing interest in quantum fluctuations and resultant 
phenomena, the Lifshitz transition has attracted rela- 
tively less attention as a quantum phase transition, in 
comparison with those of quantum critical points ob- 
served, for example in the itinerant electron magnets. 8-10 
For conventional magnetic transitions, quantum critical- 
ities emerge by reducing finite critical temperatures to 
zero, beyond which the transition disappears. On the 
other hand, as we will show in this paper, the Lifshitz 
criticality has features very different from that of the 
conventional magnetic phase transitions. For example, 
because of its topological nature, the Lifshitz transition 
does not disappear even when parameters are controlled 
beyond the point of vanishing critical temperature as we 
will clarify later. 

The Lifshitz transition occurs in noninteracting 
Fermion systems, which is characterized by the topologi- 
cal change of the Fermi surface in the Brillouin zone. Lif- 
shitz 1 classified the changes of the topology of the Fermi 
surface in two types as are illustrated in Figs. 1 (a) and 
(b). The first type is the case that the "neck" of the Fermi 
surface collapses. An appearance/disappearance of new 
split-off region of the Fermi surface is the other case. We 
call the first type of Lifshitz transition neck- collapsing 
type and the second one pocket-vanishing type. 

*E-mail:yamaji@solis.t.u-tokyo. ac.jp 
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Fig. 1. (a) Collapse of "neck" of a Fermi surface, (b) Appear- 
ance/disappearance of a new split-off region of surface. 1 



Because the topology of the Fermi surface is well de- 
fined only at T — 0, Lifshitz critical point/line, in the 
original sense for noninteracting fermions, exists only at 
T = 0. Furthermore, the transition is always continuous. 
These are clearly different from the conventional quan- 
tum critical points which originate from the suppression 
of critical temperatures as described above. The singular- 
ities seen in the Lifshitz transition are also different from 
those of the conventional quantum critical points. The 
singularities of the Lifshitz transitions can be seen as sin- 
gular parts of thermodynamic potentials, f2 s i ng . In three 
dimensions, the singular part is proportional to ICI 5 ^ 2 - 1 
The control parameter £ is the Fermi level, Ef, measured 
from the critical point value, ep] ( = Ef~ e c f . In two di- 
mensions, the singular part behaves as |£| 2 In 1/ICI f° r t ne 
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neck-collapsing type and behave as £|CI for the pocket- 
vanishing type. In fact these are derived from the well- 
known density of states in two dimensions either with the 
van Hove singularity or the band edge. Namely, the sin- 
gular part of the density of states -D s i ng behaves as ICI 1 ^ 2 
in three dimensions while in two dimensions it scales as 
In 1/| C| for the neck-collapsing type and scales as Hcavi- 
side functions for the pocket- vanishing type. Here Heav- 
iside function 9{x) is defined as follows; for x > 0, 9{x) is 
equal to 1 and, for x < 0, 9(x) is equal to 0. The singu- 
lar part of the thermodynamic potential O smg is related 
with the density of states as -D s in g = d 2 fl s i ng /d( 2 , which 
leads to the above singularity. The charge compressibility 
k and the uniform magnetic susceptibility \ are scaled 
in the same fashion as Z? s ing- It has been emphasized 
that Lifshitz transition does not follow the conventional 
Ginzburg-Landau- Wilson (GLW) scheme, 11 because the 
thermodynamic potential cannot be given by an analytic 
expansion of the control parameter £ . This violation 
of GLW scheme is in sharp contrast with the quantum 
critical point of conventional magnetic transitions. These 
singularities are listed in Table I. 
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Table I. Singularities of conventional Lifshitz transitions. The 
singular parts of thermodynamic potentials f2 s ; ng and of den- 
sities of states -D s i n g are given as functions of the Fermi level 
measured from critical points £. The charge compressibility K 
and the uniform susceptibility \ behave as D s i ng . In three dimen- 
sions, all the exponents are common between neck-collapsing and 
pocket- vanishing cases as has already been pointed by Lifshitz. 1 

Lifshitz considered the three-dimensional case and 
called the anomaly "transitions of the 2^ order" based 
on the terminology of Ehrenfest; the second derivatives 
of the electronic thermodynamic potentials have non- 
diverging square-root singularities and the third deriva- 
tives have diverging inverse-square-root singularities. 

Soon after the original remarks by Lifshitz, experimen- 
tal measurements of Lifshitz transitions were reported. 
The nonlinear pressure dependences of the superconduct- 
ing transition temperatures T c were found for thallium 
and interpreted as the effects of Lifshitz transition. 2,3 For 
Li-Mg alloys, crossover effects from the Lifshitz transition 
were also detected by thermopower measurements. 6 ' 7 

Recently, the topology of the normal state Fermi sur- 
faces of the high-temperature superconducting cuprates 
(HTSC) has attracted considerable attention. From the 
beginning of the studies on the cuprate superconductors, 
it has been well known that the Hall coefficient changes 
sign near the so-called optimum doping for the high- 
est superconducting transition temperatures. 12, 13 This is 
consistent with the topological change of the Fermi sur- 
face from a hole-like to an electron-like ones. A dramatic 
change of metallic properties from underdoped to over- 
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Fig. 2. Schematic phase diagram proposed for intermetallic com- 
pounds ZrZn2 and UGe2. 18 P and H stand for applied hydro- 
static pressures and external magnetic fields, respectively. The 
broken lines on the P-H surface stand for the "crossover" lines, 
which separates the paramagnetic states, the low moment ferro- 
magnetic phase and the high moment ferromagnetic phase (see 
text). T* means the temperature where a sudden increase of the 
magnetization occurs. 



doped regions suggests possible important role of such 
Lifshitz transitions, although its profound effect has not 
fully been elucidated. 

Among various HTSCs, the topology of the Fermi 
surface of bilayered cuprates Bi 2 Sr 2 CaCu 2 08+a (Bi- 
2212) has been intensively investigated recently by angle- 
resolved photoemission spectroscopy (ARPES). While 
some ARPES studies of Bi-2212 conclude the presence 
of large hole-like Fermi surfaces persisting for a wide 
range of dopings, other ARPES measurements suggest 
the change in the topology of the Fermi surface. Recently 
the bilayer splitting has been observed around the opti- 
mal doping. 14 The other groups also have observed the 
bilayer splitting and furthermore they have concluded 
that the change of the Fermi surface topology of the anti- 
bonding band occurs within the moderately overdoped 
region. 15 

Changes of the Fermi surface topology have also been 
proposed to be a driving mechanism of metamagnctic 
behaviors in itinerant electron ferromagnets ZrZn 2 and 
UGe2. 16 Inter-metallic compounds ZrZn 2 and UGe 2 have 
been intensively studied for decades as typical itinerant 
electron ferromagnets. Recently, these compounds have 
attracted renewed interest because of the superconduc- 
tivity and multiple quantum phase transitions. 17 The 
metamagnetic behaviors and an accompanying anoma- 
lous change of the effective mass of quasiparticles are 
observed within the ferromagnetic phase, under applied 
magnetic fields. Shown in Figure 2, discontinuous tran- 
sitions separating the phase with a relatively low mag- 
netic moment and that with a high moment are also ob- 
served. 19 One possible mechanism of such discontinuous 
metamagnetic transitions is the topological changes of 
the Fermi surfaces. If it is true, it means that Lifshitz 
transitions can become discontinuous one, which contra- 
dicts the absence of finite-temperature Lifshitz transi- 
tions predicted for noninteracting systems as described 
above. 

The electron correlations have been considered to 
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play crucial roles in various phenomena observed in the 
cuprates. In the itinerant electron ferromagnets, the elec- 
tron correlations are also important. On the other hand, 
the original idea of Lifshitz transition was based on the 
single-particle picture. The electron correlation plays no 
role in the discussion by Lifshitz. Can the electron cor- 
relations play relevant roles in the Lifshitz criticalities ? 
If the nature of the Lifshitz criticalities changes under 
the electron correlation effects, how do the criticalities 
change ? 

In the present paper, to understand electron correla- 
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Table II. Singularities of Lifshitz transitions. The singular parts 
of densities of states D s ; ng are given as functions of the Fermi 
level, C, measured from the critical points. The charge compress- 
ibility k and the uniform susceptibility \ behave as -D s i ng in 
noninteracting systems while they clearly differ at MQCP for 
the pocket-vanishing type. MQCP stands for the endpoint of 
discontinuous Lifshitz transitions. The charge compressibility n 
stands for the singularity of MQCP of the grand canonical case. 
On the other hand, the susceptibility \ stands for the canonical 
case. 

tion effects on Lifshitz transitions, we construct a phe- 
nomenological theory by assuming a mean field arising 
from a preexisting symmetry-broken order. Electronic 
interactions play roles in this phenomenological theory 
through the preexisting symmetry-broken order. We also 
perform mean-field calculations on a microscopic model 
which exhibits metamagnetic transitions. In this paper, 
we mainly consider at zero temperature, T — 0. Finite- 
temperature properties will be discussed in a separate 
publication. 

For conventional metamagnetic quantum critical 
points, the universality class has been proposed to be the 
overdamped, conserving Ising type, 20 which corresponds 
to that of the clean itinerant electron ferromagnets. This 
basically follows the conventional GLW scheme. In con- 
trast, we obtain a novel quantum criticality for the meta- 
magnetic transitions driven by Lifshitz transitions, which 
is represented by a nonanalytic free-energy expansion 
due to the singularities of the density of states. It clearly 
violates the GLW scheme. For the conventional meta- 
magnetic quantum transition, the critical exponents are 
theoretically expected to be the mean-field exponents of 
Ising model. 8-10 In contrast, we obtain different expo- 
nents at the endpoint of the first-order Lifshitz transi- 
tion at T = 0, which will be called the marginal quantum 
critical point (MQCP) in the following section, because it 
appears at the margin of the quantum critical line as we 
will see later. The difference from the conventional quan- 
tum critical endpoint within the GLW scheme arises from 
the topological origin of the Lifshitz transitions. Because 
of its topological constraint, the transition does not ter- 
minate at MQCP, but a quantum critical line at zero 
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temperature continues. Although the finite-temperature 
critical point of the Lifshitz transition may be described 
by the Ising universality class, it is deeply modified at 
MQCP. 

Results at MQCP reflect electron correlation effects. 
The charge compressibility k and the uniform magnetic 
susceptibility x, we will obtain are summarized in Tables 
II. The charge compressibility k stands for the singularity 
of MQCP of the grand canonical case. On the other hand, 
the susceptibility \ stands for the canonical case. The 
charge compressibility and the uniform magnetic suscep- 
tibility are enhanced at MQCP more strongly than the 
noninteracting case. 

In Section 2, we derive a phenomenological theory to 
elucidate the nature of the novel Lifshitz criticality at the 
endpoint of the discontinuous Lifshitz transition. Nonan- 
alytic expansions of the free energy are obtained. These 
expansions clearly violate the GLW scheme. In Section 3, 
we present numerical studies on the i-i'Hubbard model 
with an external magnetic field, which exhibits such novel 
Lifshitz criticalities. Section 4 is devoted to summary and 
discussion. 

2. Free-Energy Expansion 

2.1 Mean-field theory 

We assume that there are two bands e±(k, A) arising 
from some symmetry breaking described by order param- 
eter A. For example, A may stand for the magnetization 
within the fcrromagnctism and each band represents the 
majority and minority spin subband, respectively. The 
explicit form of these spin subbands is given as, 

e±(k,A) =£ (fc)±|A + A|, (1) 

where £o(k) is the bare band dispersion and Ao is a 
real constant, which expresses a preexisting symmetry- 
breaking field. The above band structures are obtained 
from a class of Hamiltonian H such as 

H = J2^o(k) (a\a k + b\b k ) 

k 

+ A ( a lK + b l a i) + 9^n ia n ib , (2) 

i i 

where a\(a k ) and b\{b k ) represent the cre- 
ation(annihilation) operators of an electron at wave 
number k belonging to A and B subbands, respectively. 
For the ferromagnetic systems, these subbands stand 
for the up-spin and down-spin subbands. The coupling 
constant is given by g with the number operators of z-th 
site n ia = a\a i and n ib = b\b i . The interaction term 
g J2i niafiib is rewritten as 

9 ^2 n ia n lb 

i 

kqp 

= 7v7 a l a kblb p a l b kb f p a p 

kp kp 

+ Jf s a l a q b l b k- q +p- (3) 

kqp;q^k,p 
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Then we are allowed to introduce possible order param- 
eters A , n a and n& denned by A = gN^ 1 ^2 k (a k b k }, 

n a = ^ s _1 Efe(4 a fe) and ™b = A^Efcv^fc)' where 
the average (• • • ) should be self-consistently determined, 

afterwards. With these order parameters, eq. (3) is de- 
coupled to 

g^2n ia n ib 

l a k b t b P -JfYl a l b k b l a P 

s kp 

9 £ ("*4 a fc + n^ b \ b k) ~ gN s n a n b 
k 

-J2^A*aib k + Ab{a k )+N s ^. 



9 \ - t 

kp 



In terms of the ferromagnctism, A is a transverse mag- 
netic field and Re A is the magnetization parallel to Ao- 
On the other hand, \n a — nt,\ and ImA stand for the com- 
ponents of the magnetization, which are perpendicular to 
Ao. We may take n a = nb = n/2 and ImA = , because 
the nonzero \n a — nt | and ImA are not favored in terms 
of the energy. Here n is the electron density. Then, after 
the Bogoliubov transformation, 



4 ± 



T2 alT T2 bl 



(5) 



we obtain the mean-field Hamiltonian Hmf as 

h mf = J2 ( e .( fc ' A ) + 



k,r/— ± 

N s gn 2 



A 2 

+ N S — . 

9 



(6) 



The ground-state wave function for Hmf is given as 



\* M f)= n a U n a i-\°)> w 

{k\s+(kA)<e F } {q\e-(q,A)<£F} 

where Ef is the Fermi level. Then, the mean field A is 
self-consistently given as 



9_ 

N s 



9_ 

N s 



(8) 



where (0) = (§mf\ O |$mf)- The electron density n is 
given asn = (n). 

The singularities of the density of states, which arise 
from the changes of the Fermi surface topology, are cru- 
cial in Lifshitz transitions. In two-dimensional electron 
systems, the density of states for the bare dispersion 
So(k), namely, D (s) inevitably has singularities such as a 
logarithmically diverging peak and a Heaviside step func- 
tion, which correspond to Lifshitz transitions of the neck- 
collapsing type and the pocket-vanishing type, respec- 
tively. To discuss Lifshitz transitions, we should write the 



free-energy density £ in terms of the density of states as 
(Hmf) 



- lim 

N s ^oo 



N, 



dee{D.(e,A) + D+(e,A)} 



A 2 gn 2 
9 4 



(9) 



where the densities of states for e± (k, A) are denoted by 
D±(e, A). The self-consistent equation (8) can be rewrit- 
ten as 



A= g - 



de {£»_(£, A) -D+(e, A)}. 



(4) The electron density is also rewritten as 



f 



de{D-(e, A) + D + (e, A)} . 



(10) 



(11) 



In the following derivations, we assume these band struc- 
tures (1) and the self-consistent equation (10) for simplic- 
ity. However, it captures the essential physics of Lifshitz 
transitions even for other band structures. 

We can make the Fermi level e f approach these singu- 
larities of the densities of states at £f = £p by changing 
the amplitude of the mean field A. When the Fermi sur- 
face topology changes, both the stable and metastable 
solutions of the self-consistent equation (8), for a cou- 
pling constant g — gL, are denoted by A^. The electron 
density, n, for given Al and g^ is denoted by n^. 

For small <5A = (A — A^), the free-energy density £ 
may be expanded in terms of SA, but with different co- 
efficients for positive 5A and negative 8A in general as 



aSA 



b p SA 2 



c p 5A d 



aSA + b m 5A 2 



JA 3 - 



- 0(5A 4 ) 
-0{SA 4 ) 



for SA > and 
for SA < 0. Such 



different coefficients appear because of the nonanalytic- 
ity of the density of states as we will clarify later. We 
take the same a for the both sides <5A 0, since it will 
turn out to be the same. For the canonical ensemble, the 
electron density is fixed at n = . For the grand canon- 
ical ensemble, the free-energy density £ = £ — -§§^5n 
is similarly expanded with respect to Sn = n — n^. The 
expansion by Sn gives the stability condition for the elec- 
tronic phase separation in the grand canonical ensemble. 

2.2 Canonical ensemble 
2.2.1 Neck- collapsing case 

The expansions of the free energy £ depend on the 
shape of the densities of states D±(e, A). Especially, for 
small deviation of the mean field <5A, the dispersion e+ 
shifts upwards with an amount <5A and e~ shifts down- 
wards with an amount — <5A. Then the Fermi level Ef 
also shifts to keep the density n. Below, the amount of 
the Fermi level shift is denoted by ( = Ef — sp, as is illus- 
trated in Fig. 3. Expansions of D + (e, A) and Z?_(e, A) 
around the value of Fermi level e F and the mean field 
Al determine the coefficients of free-energy expansion. 
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(a) SA=0 s 



(b) 5A>o e 




Fig. 3. (a)Schcmatic illustration of density of states D±(e, A), for 
SA = 0. (b) Schematic density of states D±(e,A), for 5 A > 0. 
Shaded areas are occupied states and the Fermi level shift, with 



an amount f = ep 
e±, namely, SA. 



ep, is induced by the shifts of dispersions 



The expansions of D + (e, A) and D_(e, A) are given as 



D + (s,A L + 5A) ~df + + df + {e- 



0{(e 



D-(e, A L + 5 A) ~d£ 



(o) 



d L 1] _{e - e L F + SA) 



+ piog In ■ 



ep + 8A\ 



0{(e 



-SA) 2 } 



(13) 



where we note that dff^ + 

j(l) 



and d^} depend on A^ while 

d K L-i ^l- an( i Piog do not f° r tnc rigid band disper- 
sions (1). Here we have used the fact that D + is analytic 

around D + (e F , Al) = df^_ and can be expanded by the 

Taylor series with the first-order coefficient d^\- On the 
other hand, £)_ has the van Hove singularity with a log- 
arithmic singularity at D_(e F , Al). The analytic part of 
Z?_ is expanded with the coefficients d^]_ and (4*1. For 
D_(s, Al), the location of the logarithmically diverging 
peak is fixed as the origin of the expansion. On the other 
hand, because of the constraint of the fixed electron den- 
sity, the coefficient of the expansion (12) depends on the 
location of the Fermi level ep and accordingly on A^ . 

When the electron density is fixed at n — til, this 
constraint determines C = e F — ep. By substituting the 
expansions of the density of states (12) and (13) into eq. 
(11), we obtain the relation between £ and S A as 



c 



-SA 1 



2d£» 



1 



Plog In 



J5A\ 



(14) 



This relation is derived in Appendix A (see eq.(A.ll)). 
From eq. (9), the expansion of £ around the Lifshitz crit- 
ical point is obtained as 

E v =const. + aSA + b n 8A 2 



SA 2 



(higher order terms) (15) 



where r\ = p,m and the coefficients a, b n and c v are given 



as 



2Ai 
11 



bp b m — 



Sg, 



9l 



-2df + + 0{5g), 



2d 



(0) J 
L + 



Plog 



(16) 



(17) 



(18) 



The higher order terms of the order 0(5g) in b p and 
b m do not play roles in the later discussion and can be 
ignored. The derivation of cq. (15) is also given in Ap- 
pendix A. The coupling constant g, which leads to Al 
is denoted by §l and 5g stands for the coupling con- 
stant measured from gv,Sg = g — gL- A striking feature 
of this expansion is that only d^i and p\ og determine the 
expansion (15). The term proportional to SA 2 / In 1 /|<5A| 
reflects the logarithmic divergence of the density of states 
(13). 21 



(12) 2.2.2 Pocket-vanishing case 



(a) SA<0 s 



(b) 5A=0 6 




Fig. 4. (a)Schcmatic illustration of densities of states D±(e,A), 
for <5A < 0. (b) Schematic densities of states D±(e, A), for SA = 
0. Shaded areas arc occupied states and the Fermi level shift, with 
an amount £ = ep — ep, is induced by the shifts of dispersions 
e±, namely, SA. 



Instead of (12) and (13), the expansions of D±(e, Al + 
5 A) for the pocket- vanishing type of Lifshitz transitions 
are given as, 



D+(e,A L + 5A) = 



9(e -ep- 5A) 



+41 (e - e 



d {0) 

a L+ 

8A) 



+ 0(s 



s Ay 



(19) 



D- (e, Al + SA) 



7(0) 



7(1) 



SA) 



+0(e-e L F + SA) . (20) 



The shifts of these densities of states induced by SA are 
illustrated in Fig. 4. Then, the condition of the fixed 
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electron density n 
( =9(-SA) 



ni, gives the following relation 



x < 



df -df + 
'df_ + df + 



5A + 2 



4mt 



(1) j(0) ■ 



- dVtd 



L- 



{dt 



4°!)- 



SA 2 + 0{SA 3 



9(6 A) (-(SA) . 



(21) 



Equation (21) is derived in Appendix A (see eq. (A. 16) 
and below) . The free-energy expansion reflects the singu- 
lar behavior of the density of states due to 9(—SA) and 
is derived in Appendix A as 



£ n =const. 



iSA 



b v 8A 2 



c 71 6A A 



where rj 



O (SA 4 ) , 

p, m and the coefficients are defined as 



(n_ L - n +L ) + 



b m = 



1 

9l 



24°>4°> 



AO) 

l L- 



2 A, 



0(5g), 



2Ai 



Sg, 



b„ = 



— + 0(Sg), 

9L 

4 4 1 M°1 3 - 



,(0) 



,(0) 



=0. 



(22) 

(23) 
(24) 
(25) 

(26) 
(27) 



Similarly to eq. (17) we ignore the higher order terms of 
0(Sg) in eqs. (24) and (25) in the later discussions. If 
c m < 0, higher order terms are needed for the stability 
of the expansion. 

2.2.3 Phase transition 

The free-energy expansions (15) and (22) have a non- 
analytic form which violates the GLW scheme. These ex- 
pansions can describe continuous and first-order tran- 
sitions depending on the sign of the quadratic coeffi- 
cient b m . For the neck-collapsing transition described 
by eq. (15), where b p = b m , the instability occurs for 
bp = b m < 0. This leads to the first-order transition 
between 5 A > and <5A < as we see in Fig. 5. On 
the other hand, the continuous transition takes place at 
a = when b p = b m > as we see in Fig. 6. For the 
pocket-vanishing transition described by eq. (22), where 
bp > b m , the instability starts first when b m becomes 
negative as is illustrated in Fig. 7. This leads to the 
first-order transition. The continuous transition occurs 
at a = if bp > b m > as is illustrated in Fig. 8. There- 
fore, for both cases, the first-order transition can appear 
only when b m < 0. Then, it is useful to introduce a di- 
mensionless factor called Lifshitz factor defined as 



1 



tl 



AO) 

1 L + 



AO) 

l L- 



(neck-collapsing) 
(pocket-vanishing) 



(28) 




Fig. 5. Free-energy density £ as a function of <5A, for the neck- 
collapsing case. For b p = b m < 0, £ has a double well structure. 
The value of 5A, namely SA P > and <5A m < give the double 
minimum of £. 



(b) a<0, b p =b m >0 





Fig. 6. (a) and (b): Free-energy density £ as a function of 8 A, for 
neck-collapsing case. For b p = b m > 0, £ has a single minimum 
at a value of 5A depending on a. For a > 0(a < 0), £ has a 
minimum at 8A < 0(<5A > 0), as is illustrated in (a) ((b)). 



b m <0 




Fig. 7. Free-energy density £ as a function of <5A for pocket- 
vanishing case. The red solid curve stands for £ p and the blue 
solid curve stands for £ m . For fe m < 0, £ has in total a double 
minimum at 8A P > and <5A m < 0. 



which is proportional to b m . The condition ?*l = is anal- 
ogous to the Stoner condition for the ferromagnetism. 
The first-order transition boundary located in the re- 
gion tl < and the continuous transition line located 
at a = 0, > meet at the point given by a = = 
at T — as shown in Fig. 9. This is nothing but the 
marginal quantum critical point (MQCP). It indeed ap- 
pears at the margin of the quantum critical line. 22 The 
appearance of the continuous Lifshitz line beyond MQCP 
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(a) a>0, b p >b m >0 



(b) a<0, b p >b m >0 



(22) that 




Fig. 8. Free-energy density £ as a function of <5A for pocket- 
vanishing case. The red solid curve stands for £ p and the blue 
solid curve stands for £ m ■ For b p > b m > 0, £ has a single mini- 
mum at a value of 8A depending on a. For a > 0(a < 0), £ has 
a minimum at <5A < 0(<5A > 0), as is illustrated in (a) ((b)). 



at T = is ascribed to the fact that the transition is a 
topological one and can not terminate. 



discontinuous Lifshitz transition line 

b<0 



marginal quantum critical point 

7 a=0,b=0 ! 

continuous Lifshitz transition line 

a=0, b >0 



o 



Fig. 9. Schematic phase diagram in plane of rj, and g. MQCP is 
the point where the conditions = and a = are satisfied. 
The thin solid line represents the continuous transition line and 
the thick solid line represents the discontinuous one. 



2.2.4 Marginal quantum criticality 

The critical exponents for MQCP can be obtained from 
the expansions (15) and (22). 

critical exponent; (3 At the first-order transition, 
levels of the two minima of £ are required to be equal. 
The jump of the mean field A at the first-order transition 
is given from eq. (15) or eq. (22) as the difference between 
the values of 6 A giving the two minima of the free-energy 
density, namely, SA p > and SA m < 0. 

In the case of the neck collapsing, this requirement is 
automatically satisfied for a = because £ is the even 
function of 5 A. Then, the first-order transition boundary 
is defined by a = and tl < 0. For the neck-collapsing 
case, the jump is given from eq. (15) as 



SA n - 5A r 



2 exp 



2d 



(0) 



fjL 



(29) 



When the exponent (3 is defined as SA p — 5A m cx I' 3 , 
it indicates (3 — > 00 (see Appendix B). On the other 
hand, for the pocket-vanishing case, we obtain from eq. 



3Ct7 



£>{rl) 



(30) 



and tl < give the condition for the two minima £ p and 
£ m being equal as Fig. 7, thus determine the first-order 
transition boundary. 

The jump for the pocket-vanishing case is given as 

12 

„3 ' 



6A % 



SA m = 



3 Cm 



9L\r L \ 
6c„ 



+ 0(ri), 



(31) 



which shows = 1. 

critical exponent; 5 At the endpoint of the first- 
order transition boundary defined by ^ = 0, (in- 
dependence of 5A determines the critical exponent 5 
from the scaling <5A cx | <5<7 1 1 Z* 5 . This relation is obtained 
from d£ p /dSA = and/or d£, n /dSA — at tl = in 
eq. (15). For the neck-collapsing case, <5A behaves as 



7(0) \.2 y \5g\' 



(32) 



24; 91 

which indicates 5=1+0 (see Appendix B). On the other 
hand, for the pocket- vanishing case, <5A is given similarly 
from eq. (22) as 



l<5AI = 




(Sg < 0) 
(Sg > 0) 



(33) 



which shows S = 2 for Sg < and S — 1 for Sg > 0. 

critical exponent; 7 Around the endpoint = 0, 
the uniform susceptibility x = (d 2 £ /dSA 2 ^) behaves 
as 

X = g L rl\ (34) 



for both of the neck-collapsing and the pocket-vanishing 
cases. The relation x ^ r L 7 with 7 = 1 is satisfied. 
However, for the the pocket-vanishing case, in the side of 
SA > 0, the uniform susceptibility x becomes a constant, 
gL, resulting in 7 = 0. 

Next we consider the scaling of x as a function of SA, 
strictly at = 0. Near MQCP of the canonical ensem- 
ble, the susceptibility x scales, as SA approaches 0, as 

■ In , „ „ , (neck-collapsing) 



X 



SA\ 



1 



&c m \SA\ 



(pocket-vanishing, Sg < 0) 



(35) 



The critical exponents for both of the neck-collapsing 
and the pocket-vanishing cases are shown in Table III. 
Each set of exponents satisfy a scaling relation 7 = (3(6 — 
1) (see Appendix A). We note that the universality class 
for the pocket-vanishing type is the same as that of the 
metal-insulator transition. 22 ' 23 

2.3 Grand canonical ensemble 

For the grand canonical case, the instability towards 
the electronic phase separation exists. The endpoint of 
the electronic phase separation is nothing but MQCP in 
the grand canonical ensemble. We clarify the instability 
by calculating the free-energy expansion with respect to 
6n, namely £. For the coupling constant fixed at g = gL, 
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NC 


PV 

(Sg < 0) 


PV 

(5g > 0) 


p 


oo 


1 


1 


7 


1 


1 





5 


1 + 


2 


1 



Table III. Critical exponents of MQCP. NC and PV stand for the 
neck-collapsing and pocket-vanishing types, respectively. 



we introduce an amount of the chemical potential shift, 
Sfj, as a control parameter as 



£ = £ 



(36) 



Then the uniform charge compressibility k, which is 
given as 

-l 

(37) 



d 2 £ 



dSn 2 



Sn=0 



exhibits a significant feature at MQCP as shown below. 

2.3.1 Neck- collapsing case 

From eq. (C.15) in Appendix C, the electron density 
(C.l) gives the relation between £ and Sn as 



c 



2s 2 



i 



1 

2rz, V 9LP\o S r L In w 



(38) 



where we assume ^ 0. Around MQCP of the grand 
canonical ensemble, the free-energy density £ can be ex- 
panded with respect to small Sn as is mentioned in §2.1. 
The free-energy expansion of £ is given as 



f — 



£,, - ^-Sn - Sp:Sn 



"0 2 

const. + aSn -+ 



b v Sn 2 



Sn" 



where n = p,m and the coefficients are defined as 
a = —6fi, 

9l fl-r L 



bra — bp 



= c P = - (4- + 4°| + Piog - 9 L *) 



2pi og r 2 L ' 



(39) 

(40) 
(41) 

(42) 



The derivation of eq. (39) is given in Appendix C (see 
cq. (C.17)). 

2.3.2 Pocket-vanishing case 

From eq. (C.25) in Appendix C, the electron density 
(C.l) gives the relation between £ and Sn as 



to make the pocket Fermi surface emerge always for Sn < 
0. The free-energy densities £ p for the positive Sn and 
£ m for the negative Sn may be separately expanded with 
respect to Sn as 

£ v =£ v — L Sn — 5/iSn 

+aSn + b v 8n + c n 5n \ (45) 
where n — p,m and the coefficients are defined as 
a = —5/j,, 



b m = 



1 



9L 



i J.0) AO) 

qT ~ a L- a L+ 



24°) 



. d (0) 1 2d (°) d (°) ■ 



1 1 



(46) 



(47) 



6 (4°1+4°1) 3 



x { (2r+ 3 + 3r-r+ 2 ) + (ir^ + 3r+r~ 2 ) d£j.} 



5 1 



(0) \ 



(48) 



where 5fx = Sfi ■ sign(4°l — 4°+), r + = g L x — c?^°]_ and 
r~ = g^ 1 — d^_. The derivation of eq. (45) is given in 
Appendix C. 

2.3.3 Phase transition 

For the neck-collapsing case, the charge compressibil- 
ity k is scaled as 

«=--^, (49) 

along the Lifshitz critical line. The derivation of eq. (49) 
is given in Appendix C (see eq. (C.20)). Here, the in- 
equality tl < 1 is always satisfied, leading to k < 0. The 
negative value of the charge compressibility k around the 
continuous transitions for > 0, means the instability 
toward the electronic phase separation. Then, the charge 
compressibility k is negative along the Lifshitz critical 
line, except for the case g^ = 0. When approaches 
zero, from cq. (28), we see that tl — * 1, resulting in 
k — ► oo from eq. (49). In fact, the phase separation 
terminates just at g^ = 0, which means that MQCP is 
located at <?l — 0. At g^ — 0, the charge compressibility 
behaves as 



1 

9L 



9L 



7 (0) 



Sn. 



(43) 



Around MQCP of the grand canonical ensemble, the free- 
energy density £ can be expanded with respect to small 
Sn as is mentioned in §2.1. For the sake of unified treat- 
ment, we introduce a modified form of Sn as Sn defined 

by 

Sn — Sn ■ sign(c4°l — 



(44) 



(50) 



as is derived in Appendix C (see eq. (C.21)). The free- 
energy expansion of £ around MQCP, namely around 
g = is given as 

1 Sn 2 



£ ~const. — SfiSn - 



~const. — (tyipiogCln 



2piog In Wl 

1 , Plog 



1(1 



C 2 ln 



icr 



(51) 
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which exhibits the same singularity mentioned in Table 
I. Here we use the relation between Sn and ( at gL = 0, 
namely Sn ~ pi og £lnl/|£|, which is obtained from eq. 
(38). Eventually, the marginal quantum criticality turns 
out to be nothing but the criticality of the noninteracting 
case. The endpoint of the phase separation only exists at 
gL = 0, as is illustrated in Fig. 10 (a). 

For the pocket-vanishing case, for SA = 0, k is given 

as 



2dV 



--9L- 



91 



7(0) 



(52) 



For the derivation of eq. (52), see Appendix C (see eq. 
(C.32)). Here the inequality r L > 1 - g L \J d^ldfl 

always holds because the geometric average \J d^_df+ 

is greater than the harmonic average 2d^_d^]_ / (d^l_ + 

df+). Then the electronic phase separation inevitably 
occurs around MQCP of the canonical ensemble, where 
= 0. The endpoint of the electronic phase separation, 

which is defined by 1 — gL \J d^_d^\_ = 0, also exists on 
the continuous Lifshitz transition line. The endpoint of 
the electronic phase separation, which is located at a = 
and b m = 0, is MQCP in the grand canonical ensemble. 
The Lifshitz factor for the grand canonical ensemble, r£, 
is given as 



= l-g L y/^ 0) 



(53) 



which is proportional to b m . The instability towards 
the electronic phase separation signaled by the nega- 
tive charge compressibility exists around MQCP of the 
canonical ensemble. We should note that the instability 
towards the phase separation only exists at the one side 
of the Lifshitz transition boundary of the canonical en- 
semble, where the pockets of the Fermi surface exist, as 
is illustrated in Fig. 10 (b). 

2.3.4 Marginal quantum criticality 

The critical exponents 0, S and 7 for MQCP of the 
grand canonical ensemble are obtained as follows. 

critical exponent; [3 The jump of the electron den- 
sity n at the first-order transition is given as the differ- 
ence between the values of Sn at the two minima of the 
free-energy density, namely, 5n + > and Sn- < 0. The 
exponent (3 is obtained from the scaling 5n + — <5n_ cx 
|frm|' 3 oc IrX] 13 . For the pocket-vanishing case, we obtain 
from eq. (45) that 



a = — ■ 



4c„ 



+ 0(b m ) 



(54) 



and b m < 0, namely tl < determine the first-order 
transition boundary. The jump of the electron density is 



(a) neck-col lapsing case 



(b) pocket-vanishing case 




9 CM0CP 



jaocp _ 

a =0 b, n =t^0 



GMOCP 

a =0, b m = 




Fig. 10. (a)Schcmatic phase diagram for neck-collapsing case. 
(b)Schcmatic phase diagram for pocket-vanishing case. In ad- 
dition to MQCP of the canonical ensemble (CMQCP), MQCP 
of the grand canonical ensemble (GMQCP), defined by a = 
and b m = 0, appears. Within the shaded areas, the charge com- 
pressibility k becomes negative and the instability towards the 
electronic phase separation exists. For the pocket- vanishing case, 
the region where k becomes negative, only exists at the one side 
of the Lifshitz transition boundary, defined by a = 0. 



given as 



Sn + — Sn_ ~ 



8c rs 



8c m |6 p | 



(55) 



which indicates (3=1. 

critical exponent; 8 At the endpoint of the first- 
order transition boundary defined by b m = 0, Sn- 
dependence of Sn determines the critical exponent S from 



the scaling Sn cx 



Sn 



1/6 



This relation is obtained from 



d£ m /d8n = and d£ p /d8n — at b m = r~L = in eq. 
(45). For the pocket-vanishing case, the behavior of Sn 
is given as 



\8n\ 




{Sn > 0) 



(56) 



which shows 5 = 2 for Sfi < and S = 1 for 5u > 0. 
critical exponent; 7 Upon approaching the endpoint 
b m = r~~L = 0, the uniform charge compressibility k = 

(^d 2 £/dSn 2 ^J is given from eq. (45) 



as 



(2b,) \ 



For Sn < 0, the relation k 



(57) 



with 7 = 1 is satis- 



fied. On the other hand, for Sn > 0, k becomes constant 
value d L _ resulting in 7 = 0. In addition, as a function 
of Sn the charge compressibility n diverges as 

1 



Qc m \Sn\ ' 



(58) 



strictly at b m = for Sn — > 0—. The scaling (58) is 
equivalent to Sn oc /i 1 /" 5 with the chemical potential n = 
d£/dSn and the exponent S = 2. In fact the exponent S 
derived in eq. (56) should be the same as S derived here 
because a = —S/j, is nothing but the chemical potential 
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These critical exponents for the pocket-vanishing case 
of the grand canonical ensemble are the same as those 
of MQCP of the canonical ensemble, which are shown in 
Tabic III. 

3. Numerical Studies 

In the present section, we demonstrate numerically the 
marginal quantum criticalities of the canonical ensem- 
ble. Here, we consider a simple model for ferromagnets, 
namely t-t' Hubbard model on square lattice defined by; 

H =^Te a (k)cl a c ka + Uj2n lA n hl , (59) 

ka i 

So(k) = — 2t (cos k x + cos k y ) + At' cos k x cos k y , (60) 

where c\ a (c k(J ) represents the creation(annihilation) op- 
erator of an electron at wave number k with spin a. The 
nearest-neighbor and next-nearest-neighbor transfers are 
denoted by t and t', respectively. The onsite interaction 
is given by U and the number operator of the i-th site 
is n ia — c\ a c i(r . It has been suggested, by a projection 
quantum Monte Carlo study, that, for It' jt > 0.55, the 
t-t' Hubbard model exhibits weak coupling ferromagnetic 
order, instead of anti-ferromagnetic one at the filling cor- 
responding to the van Hove singularity. 24 In addition to 
the t-t' Hubbard Hamiltonian (59), we introduce a trans- 
verse magnetic field H± as follows; 

^E^Va+^W), (61) 

i 

where the Lande's g-factor and the Bohr magneton fis 
arc included in the notation of H±. Under a nonzero 
magnetic field H±, the magnetization 

m^^N-^icl^) (62) 
k 

always remains finite. 

To study the mctamagnetic transitions within the 
present model, we use the mean-field approximation, and 
the terms of the on-site Coulomb repulsions are decou- 
pled with the mean field m± as follows; 

kq 

- - E ( Um ± c U c ki + Um i c li c ki) 
k 

+N s U\mj_\ 2 . (63) 

Then, we obtain the quasiparticle dispersions of the spin- 
polarized states e±{k) as, 

e±(k)=e (k)±\Um± + H ± \, (64) 

after the Bogoliubov transformation, 

«L = 7f C feT T ^f4i- (65) 
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Here we assume that m± ■ H± > 0. The electron density 
n and m± are calculated by the following conditions; 

n =]v-E{/( e -( fc )) + /( e +( fc ))}> ( 66 ) 

s k 

m^=^J2{f(e.(k))-f(e + (k))}, (67) 

s k 

where f(x) = (1 + exp [(x — ^i)/T]) _1 and the chemical 
potential fx and the magnetization mj_ must be deter- 
mined self-consistently for the given electron density n. 

In analyzing the t-t' Hubbard model with the exter- 
nal magnetic field by the mean-field approximation, we 
take t' = OAt as a typical value and leave n, U and H± 
as control parameters. To avoid the instabilities toward 
commensurate antiferromagnetic orders, at least within 
the framework of the mean-field theory, we consider the 
low-electron-density region n < 0.6. The overall phase 
diagram is shown in Fig. 11. For values of the electron 
density n > 0.29, the Fermi level is located below the van 
Hove singularity of the majority spin sub-band £-(k) at 
U = 0. As U increases, the Fermi level approaches the 
van Hove singularity and a neck-collapsing Lifshitz tran- 
sition occurs. Above the neck-collapsing Lifshitz critical 
value of U, we further observe the pocket-vanishing Lif- 
shitz transition. 

On the other hand, for the lower electron density n < 
0.29, the Fermi level is always located below the van Hove 
singularity of the majority spin sub-band e_ (fc) and only 
pocket-vanishing Lifshitz transitions can occur. These 
pocket-vanishing Lifshitz transitions at lower electron- 
densities are nothing but the transitions between par- 
tial and complete ferromagnetic states. For both of the 
neck-collapsing and the pocket-vanishing cases, the Lif- 
shitz transitions contain both discontinuous and continu- 
ous parts. MQCPs appear at the boundaries between the 
first-order and continuous transition lines, as are shown 
in Fig. 11 as open circles. 

Let us consider at the electron density n = 0.4, for 
example. Then the neck-collapsing Lifshitz transitions 
occur at the boundary between (c) and (d) in Fig. 11. 
It becomes discontinuous for relatively small values of 
H±. The jump of the mean field at the first-order tran- 
sition rapidly decreases to zero with increasing H±, as 
is shown in Fig. 12. The concave behavior of the jump 
(5A + — <5A„ as a function of H±/t indicates that the crit- 
ical exponent (3 should be larger than 1. As shown in Fig. 

13, near MQCP, <5A behaves as 8 A oc 8glnl/\5g\. This 
scaling is consistent with eqs. (29) and (32) indicating 
(i = oo and (5=1 + 0. 

On the other hand, the density dependence of 5A + — 
<5A_ for the pocket-vanishing transition is shown in Fig. 

14. Since the electron density measured from the critical 
electron density at MQCP Umqcp is scaled linearly with 
tl, it leads to (3 — 1. As is shown in Fig. 15, the critical 
exponent for the pocket-vanishing case is given by 5 = 2. 
These exponents are consistent with eqs. (31) and (33). 

4. Summary and Discussion 

In the present paper, electron correlation effects on 
Lifshitz transitions have been investigated by assuming 
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Fig. 11. Upper panel, (a),(b),(c) and (d): Shapes of Fermi sur- 
face corresponding to the regions in the phase diagram in the 
lower panel. Lower panel: Phase diagram of Lifshitz transitions 
at t' /t = 0.4 and H± = 0.06i. The thin black curves repre- 
sent the continuous transition lines of the pocket- vanishing case. 
The thin blue curve represents the continuous transition line 
of the neck-collapsing case. The thick red curve represents the 
first-order transition line. The vertical arrows indicate MQCPs. 
The phase boundary between (a) and (c) as well as between (b) 
and (d) represents the pocket-vanishing transition. On the other 
hand the boundary between (c) and (d) as well as (a) and (b) 
represents the neck-collapsing transition. 



the presence of a broken-symmetry order. Through the 
interplay with the order parameter amplitude, electron 
correlation causes the first-order transition in contrast 
with the noninteracting case. The first-order transition 
boundary appears in a part of the parameter space while 
a continuous transition boundary also exists. These two 
boundaries meet at the marginal quantum critical point 
(MQCP) as seen in Fig. 9. The continuous boundary at 
T = appears because this transition is a topological 
one. The topological nature deeply influences the quan- 
tum criticality and the criticality does not follow the 
Ginzburg-Landau- Wilson (GLW) scheme. 

Let us summarize the parameters to be controlled here. 
The Fermi level £ measured from the critical point is 
the control parameter in the noninteracting case. In the 
present interacting system, the primary control parame- 
ter may alternatively be taken as <5A, the order parame- 
ter of the broken symmetry. Then £ is determined from 
5 A self-consistently. In terms of the experimental con- 
trol, SA is self-consistently determined from an external 
control parameter such as the interaction parameter Sg 




0.08 



Fig. 12. Jump of mean field, <5A+ — <5A_, as a function of H±/t 
at fixed electron-density n = 0.4. As H±/t increases, the jump 
rapidly decreases to zero as a concave function of Hj_/t. 




log I Vsg I 



Fig. 13. Ratio of 5A to 8g plotted as a function of lnl/|<5g| near 
MQCP at U ~ 1.478t,n ~ 0.463 for H± = 0.064. The linear 
behavior of <5A indicates that <5A <x 8g In 1/|<S<?| is satisfied. 




0.005 0.010 015 0.020 0.025 

Fig. 14. Jump of mean field, 5A4. — <5A_, as a function of \n — 
"mqcpI a t H± = 0.064, where the critical electron density at 
MQCP is denoted by rtMQCP ■ 



measured from the critical point, which may be directly 
expressed by pressure, for example. Another controllable 
parameter is the Lifshitz factor as is seen in Fig. 9. 
The Lifshitz factor tl is considered to be controlled by 
the external magnetic field H± and/or the electron den- 
sity n. 

We first summarize the case of the canonical ensemble. 
In the presence of interactions, the singularities of the 
continuous Lifshitz transitions seen in the susceptibilities 
X are different from the noninteracting case shown in 
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Fig. 15. Log-log plot of <5A| vs \8g\ around MQCP of pocket- 
vanishing case for negative Sg and H± = 0.06t. The relation 
|<5A| oc l^gl 1 / 2 is obtained. 



Table I. For r L 
f In 



X <x 



| SA I 

|<5A|-' 



0, the susceptibilities x behave as 

(neck-collapsing) 
(pocket-vanishing) 



(68) 



which corresponds to eq. (35). In terms of the control 
parameter Sg, the susceptibility \ behaves as 



X oc 



In J- 

Sg 

Sg' 1 ' 2 



(neck-collapsing) 
(pocket-vanishing) 

In terms of tl, the susceptibility \ behaves as 

„-i 



X 



(69) 



(70) 



In contrast to the noninteracting case, where the suscep- 
tibilities do not diverge, the marginal quantum critical- 
ity exhibits the diverging fluctuation, especially for the 
pocket-vanishing case. We also note that the shift of the 
Fermi level £ behaves as 



-8A (neck-collapsing) 



CCK 



d[ 0) -d 



1.0) 



— } (of 8 A (pocket-vanishing, 6 A < 0) , (71) 

—SA (pocket-vanishing, SA > 0) 

when 8 A — ► 0. 

Next we summarize the grand canonical case. For 
the neck-collapsing case, the electronic phase separa- 
tion always occurs around the Lifshitz critical line, in 
the presence of the interaction. On the other hand, for 
the pocket- vanishing case, the endpoint of the electronic 

phase separation, which is defined by g \J d d — 1, ex- 
ists. Then the charge compressibility n diverges as 

k oc \5A\- 1 oc \8n\~ 1 cx \5(j,\~ 1/2 . (72) 

We also control the divergence of k around grand canoni- 
cal MQCP by changing coupling constant Sg, which turns 
out to be proportional to Sji. In terms of rj,, k is scaled 
as 



K OC 



(73) 



where = 1 — g y d a d a . At the grand canonical MQCP 

where f£ — ► 0, k diverges as r^ -1 . In this case, the shift 
of the Fermi level ( behaves as 

C oc Sn, (74) 

when Sn — ► 0. 



The macroscopic electronic phase separation is pre- 
vented by the charge neutrality, in the metallic com- 
pounds. However, the instability towards the phase sep- 
aration itself exist intrinsically and the uniform state is 
unstable anyhow. This instability around MQCP of the 
grand canonical case competes with the charge neutral- 
ity. Then the competition would result in the electronic 
inhomogencity. 25 The existence of the electronic inhomo- 
geneity has been discussed for the HTSCs 26 and the man- 
ganites 27 and has attracted increasing interest. The elec- 
tronic phase separation induced by the two-dimensional 
van Hove singularity, which is nothing but the phase sep- 
aration induced by the neck-collapsing Lifshitz criticality, 
was also discussed in the context of the multiple metam- 
agnetic transition observed in the bilayer Ruthenates. 28 

For both of the canonical and grand canonical ensem- 
bles, we obtain the critical exponents of MQCP, which 
are different from those of the mean-field Ising universal- 
ity, as is summarized in Table III. It is necessary to fur- 
ther examine the validity of these exponents, since these 
critical exponents are obtained by the mean-field theory, 
which ignores spatial and temporal fluctuations around 
critical points. For the neck-collapsing case, because of 
the logarithmic singularity, it may be difficult to deter- 
mine the critical exponents of MQCP in real systems 
although the first-order jump may be observed. On the 
other hand, for the pocket-vanishing case, the power-law 
singularities may easily be observed. 

The free-energy expansion and the mean-field expo- 
nents for MQCP of the pocket-vanishing case may be 
justified beyond the mean-field approximation. This is 
because the dynamical exponent z is expected to be 
equal to 4. In two dimensions, the effective dimension of 
the marginal quantum critical phenomena of the pocket- 
vanishing case, d + z, may be at the upper critical di- 
mension of the y> 3 -theory, namely d c — 6. 23 The details 
will be reported in a separate publication. 

Although we consider the ferromagnetic band disper- 
sion (1) and the corresponding self-consistent equation 
(10) for simplicity, the essential physics of the marginal 
quantum criticality is retained if individual band disper- 
sions are considered in other possible cases. The ferro- 
magnetic model treated in this paper can in fact also be 
mapped to another model, namely t-t' Hubbard model on 
two layers which contains homogeneously coupled two- 
dimensional square lattices. Then the spin degrees of 
freedom in the ferromagnetic model are mapped to the 
layer degrees of freedom. These two degrees of freedom 
are referred to as A and B subbands in §2. The exter- 
nal magnetic field is mapped to an inter-layer transfer 
integral. Additional spin degrees of freedom exist in the 
bilayer model. When the magnetic orders are absent, the 
susceptibility x here represents the inter-layer fluctuation 
of the electron density. This bilayer model is expected 
to have common features to the layered two-dimensional 
metals like Bi-2212. 
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Appendix A: Derivation of free-energy expansion 
of canonical ensemble 

neck-collapsing case For the band structure 
e±(k,A) = £o(k) ± (A + A), the density of states 
D±(e 1 Al + SA) can be replaced by 

D ± (e-e L F TSA) = D ± {e,A L + SA), (A.l) 

where the right hand side is defined in eqs. (12) and (13). 
Here, we assume A • A > 0. Then the condition for the 
fixed electron density (11) can be rewritten as 

/C+«5A f-C-SA 
dC-D_((_) + d( + D + ((+), (A.2) 

-oo J — oo 

where (± is given as C± = £ — £p T SA. The critical 
electron density at the Lifshitz transition tll is defined 

as 



= [ dC-S-(C-)+ f d( + D + (( + ). (A. 

J — oo J — OO 



3) 



Then the condition n — results in 

=n — riL 

rC+SA r(SA 

= / dC-I>-(C-)+/ d( + D + (C+). (A.4) 
Jo Jo 

For the neck-collapsing case, the expansions of the den- 
sity of states (12) and (13) are rewritten as 

D + (C+) =4°1(Al) + 4 1 1(Al)C+ + O (Cl) , (A.5) 
D_(C_) =4°1 + Plog In + 4°lC- + O (C_) • (A.6) 

By substituting the expansions (A.5) and (A.6) into 
(A.4), we obtain an approximate formula for the con- 
straint of the fixed electron density as 

Plos ln ](T5A\ (C + SA) + K" + Aog ) (C + SA) 

+d { °l((-SA) + 0(5A 2 ) =0, 

(A.7) 

up to the order of <5A 2 . To obtain the solution for eq. 
(A.7), we may take the following form of £, 



C = -JA(l-ft), 



(A.8) 



where TZ is a function of SA and is defined by this rela- 
tion. From eq. (A.7), it is apparent that SA is propor- 
tional to £ except for a logarithmic correction so that TZ 
may have an amplitude of order unity. Now cq.(A.7) can 
be rewritten by using (A.8) as 



1 



p log S Aft In — + (4°{ + Plog )SATZ 

-d ( "l(28A- 8AU) + 0(5A 2 ) =0. (A.9) 



We can drop linear terms of 5A1Z to discuss the 
asymptotic behavior, since 5A1Z is higher order than 
SAlZluj^-. Then eq. (A.9) is simplified as follows; 



TZ 



2d 



(0) 

L+ 



1 



Plog In 5An 



(A.10) 



The asymptotic solution of eq. (A.7) is accordingly ob- 
tained by solving (A. 10) iteratively as 



c 



-5A 1 - L+ 



(A.ll) 



which completes the derivation of eq. (14). The presence 
of the logarithmic divergence due to the van Hove singu- 
larities makes the asymptotic behavior determined only 
by 4+ and pi og . The free-energy density at T = given 
by eq. (9) is rewritten as 

K+SA 

£o= dC_(C--<fA)£>_(C-) 



+ 



d( + (( + + 5A)D + (( + )+ [ -^- 

oo 9 



5A) Z 



-E-L + E+L 

K+sa 



(n-L - n +L ) 5A + 
dC- (C--«A)£>_(C_) 



(A L + «5A) Z 



o 



+ 



C-«5A 



dC+(C+ + 5A)^ + (c+), 



(A.12) 



where the constant terms E±l and n±L are defined as 
£ ±L d( ± C±D ± (C±) and n ±L eee d(±D ± ((±). 

By substituting the expansions (A.5) and (A.6) and 
the relation (A.ll) into (A.12), the expansion of £o 
around the Lifshitz critical point is obtained as eq. (15). 

pocket-vanishing case For the pocket-vanishing 
case, the expansions of the density of states (19) and 
(20) are rewritten as 



D + (C+) =0(C+) 



d_(c-)=4-+<*l-(c-) + 0(c-) 



4°l+4j(c + , 



L+ 
(1) 



+ 0((l) , (A.13) 



(A.14) 



"L- 1 L- 

By substituting the expansions (A.13) and (A.14) into 
(A.4), we obtain an approximate formula for the con- 
straint of the fixed electron density as 

o=(4°l + 4 l)c + (4 l-4>A 



4|) (C 2 + SA 2 ) 



(41 



d ( ll)(SA + 0(5A 3 ), 

where we assume £ — SA > 0, otherwise the constraint of 
the fixed electron density is given as 

h d L-(t + SA ) 2 - From °q- ( A -!5), the relation 



(A.15) 



4 0) -4°l 4 1) 4°l 



41 + 4+ 



7(0) 



+ O (SA 3 ) 



(4°1+4°1) 1 



^+^°- SA 2 



(A.16) 
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—SA > 0, otherwise £ = — <5A is 



for (-SA = 
obtained. 

By substituting the expansions (A. 13) and (A. 14) and 
the relation (A. 16) into (A. 12), the expansion of £o 
around the Lifshitz critical point is obtained as eq. (22). 

Appendix B: Critical exponents of neck- 
collapsing case 

Because of the van Hove singularity, the critical ex- 
ponents are not well-defined on the marginal quantum 
critical point of the neck-collapsing case. The relation, 
Sg cx tJA/lni, has to correspond to Sg cx SA S . Here, 
to focus on the slow convergence of 1/ In as 5 A ap- 
proaches zero, we use a symbolic notation as 1/ In cx 
SA +a . This notation represents that 1/ln^ converges 
to zero more slowly than 5A P , where the exponent p takes 
the limit +0. From eq.(29), SA depends on < as the 
following essentially singular form, 

\SA\ cx cxp [-const, x Ir^p 1 ] , (B.l) 

when tl approaches zero along the first-order transition 
boundary. Equation (B.l) yields the following relation, 

'^hrW < B " 2 ) 

Then, by using the above notation, the relation \tl\ oc 
|c5A| +0 holds symbolically. It may be rewritten as |<5A| cx 



1/+0 



symbolically. On the other hand, because the 
relation (lnl/^A)" 1 cx 5g/6A or the symbolic form 
|<5A| +0 cx Sg/5A holds from eq. (32), another relation 
is obtained as 



ItlI oc 



6 A' 



(B.3) 



which corresponds to a contract +0/ + = 1 for the sym- 
bolic notation, because of the symbolic relation \5A\ 



+o 



oc 



\\r L \ 1/+Q \ cx dg/SA. Then critical exponents are ex- 
pressed symbolically as 



0=-^, 5 = 1 + 0. 



(B.4) 



By noting the fact that the critical exponent 7 is equal 
to 1, the scaling relation, (3(6 — 1) = 7 holds with the 
contract +0/ + 0=1. 

Appendix C: Derivation of free-energy expansion 
of grand canonical ensemble 

To calculate the free-energy expansion for the grand 
canonical ensemble, we assume the difference between n 
and til a,s Sn = n — n^. From eq. (A. 4), Sn is given as 

r(+8A fC-SA 

Sn = dC-D_(C-)+ d( + D + (( + ), (C.l) 

Jo Jo 

where S A and £ are functions of Sn. For the fixed cou- 
pling constant g = gL, the mean field 5 A satisfies the 
following self-consistent equation 



A L + SA = 9 4 



dC-D-(C-) 



f 



-SA 
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where A^ is given as 



A, = 



9l 



{ f dC-^-(C-) - [ dC+D + ((+)} ■ (C.3) 

w— 00 J — CO ) 



Then SA is given as 
SA = 



f y^ +SA d(_D_(C-) £ dA d( + D + (C + )^ .(C.4) 

neck-collapsing case The free-energy density is 
given as 

K+SA 

£= d(- (C- - SA) D_(C-) 

Jo 

rCSA 

+ dC + (C++SA)D + (C+) 
Jo 



SA gn 2 
9 4 



+ ^(C 2 - SA 2 ) In -i— + £^(C + SA)(( 3SA) 



\C + SA\ 4 



+ 6 A + 9± +0{ 5A*). 
9 4 



(C.5) 



By substituting the expansions (A. 5) and (A. 6) into 
(C.l) and (C.4), we obtain the following relations 

Sn = (df_ + Ploe ) (C + SA) + Plog (C + SA) In 



+ d^ + ((-5A) + 0(SA 2 ), 



\C + SA\ 
(C.6) 



9L 



SA =(df_ + Ploe ) (C + SA) + Ploe (C + SA) In 



1 



- • -(0) 

d ( £>((-5A) + 0(5A' z ). (C.7) 



|C + <5A| 



From eqs. (C.6) and (C.7), we obtain the relation which 
does not include logarithmic terms as 

Sn--5A = 2d<£l((-SA)+0(5A 2 ). (C.8) 
By solving (C.8), SA is given as 

SA 



Si L + 

Then we can eliminate SA from (C.6) as 



(C.9) 



Sn 



, ,(0) . ^ g L Sn + 2r L ( 
(d> L !_ + piog) — 

gSn + 2r L ( 2s 
+Pl ° S 2s lU \g L Sn + 2r L (\ 
+d W Z 9jJn + 2C +0{§n% (C1Q) 



2s 




where = l — 2gi,d^ and s = l — g^d^. To obtain the 
solution for eq.(C.lO), we may take the following form of 



Sn = -^(1+S) 
9l 



(C.ll) 
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where S is a function of £ and is defined by this relation. 
The relation (C.10) can be rewritten by using (C.ll) as 



2?-lC „,(o)> r L (S s 
2d L+( - Piog In ■ 



(C.12) 



9 " ' ' ■ s \r L (S\ ' 

where the terms, which are order of (S, are dropped. 
Then the asymptotic form of S is given as 

2 s 2 1 



9LP\ og r L In 1 



ICI 

Then, eq. (C.ll) can be rewritten as 



Sn ~ — 



2r L ( 

9L 



1 + 



9LP\ og r L In 1 



ICI 



We can also solve eq. (C.14) iteratively as 



c 



9lSu 
' 2r L 



1 



2s 2 



9LP\o e r L In w 



(C.13) 



(C.14) 



(C.15) 



where we assume ^ 0. By substituting (C.15) into 
(C.9), we obtain SA as a function of Sn as 



SA 



9L_ Sn 9Ld { l ] + s Sn 



2r L p\o % r 2 L ln^ 



(C.16) 



By substituting (C.15) and (C.16) into (C.5), we obtain 
the free-energy expansion with respect to Sn as 

_ Sn 2 



f - 



9L n L 



Sn ~ const. + b v 5n 



'In 



(C.17) 



where n = p,m and the coefficients are defined as 



brn — bp — 



9L ( 1 - r L 



(C.18) 



S9L 



C m = C p = - (4°1 + 4+ + Aog -9L 1 ) g 2 ' 

l0g (C.19) 

From eq. (C.17), the uniform charge compressibility k is 
given as 

_i d 2 



dSn 2 



,g g 
"2 2r L 



2s 2 



5Piogr-L In w 



(C.20) 



In this case, from the inequality = 1 — <7l4+ < 1, 
the charge compressibility k is negative along the Lifshitz 
critical line, except for the noninteracting case. For 5 = 0, 
the charge compressibility behaves as 

« = Piogln — ^. (C21) 



pocket-vanishing case For £ — SA > 0, the free- 
energy density is given as 

r(+6A 

£= dC-(C--«A)D_(C_) 
Jo 

+ / dC+(C+ + ^A)^ + (c+) 

Jo 

(5A an 2 
+ — + V' 



+ 



-i(41-41)^ + ^ + -^ + 0( ^). 

(C.22) 

By substituting the expansions (A. 13) and (A. 14) into 
(C.l) and (C.4), we obtain the following relations 



Sn =41 (C + SA) + ^-(( + SAf 



7(1) 



+ 4°| (C - <5A) + -|±(C - 5 A) 2 + 0(SA 



7(1) 
L + , 



(C.23) 



3L 



Ai) 

SA =4°1(C + <5A) + -^(C + <5A) 2 



7(1) 



- 4| (C - <5A) - -|±(C - SA) 2 + 0(SA 3 ). 

(C.24) 

By solving eqs. (C.23) and (C.24), we obtain ( and SA 
as functions of Sn as 



df-+df + r n 2^(0) +d (0)y 



(C.25) 



<5A: 



lrf 



(0) ,(0) 



L " L ~ UL + i S n+- - 

2 4°1+4V 2 (41 + 41) 



CA -Sn 2 



3 r 3 ' 



(C.26) 



2d (o) d (o) 

where the coefficients are defined asrE- L ~ L+ 



9l 

C„ = 41 (£ - 41) 3 + 41 (£ - 41) 3 aud C A = 

-dj^ - 41) 2 + 414°1 {-k - 41) 2 . Here, 
we assume r ^ 0. 

By substituting (C.25) and (C.26) intoJC.22), we ob- 
tain the expansion of £ v with respect to Sn defined in eq. 

(44) as 

£ n — 9 L ™ L Sn ~ const. + b rj 5n + c ri 5n , (C.27) 
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where r\ = p,m and the coefficients arc dcfinecl as 



d L- d L+ 



2 w<°) 



d^>_ + d 



(0) 



2J 



(0) ,,(0) ' 



1 1 

w 

1 



(C.28) 
(C.29) 



9) 
10) 
11) 

12) 



13) 
14) 



6 (4°! + 4°i) 3 r3 



+ 2 



7(0) 



9L 



5 1 



+ 3 



6^0)3 



1 3 

1 - l9L d L- 




Then the uniform charge compressibility k around Lif- 
shitz critical points is given as 

d 2 L 



On 2 

n—n^+Sn 

~ 2b n + &c v 8n. 



(C.32) 
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